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Abstract. The recent observation by the PAMELA satellite of a rising positron fraction up to ~ 100 GeV has triggered 
a considerable amount of putative interpretations in terms of dark matter (DM) annihilation or decay. Here, we make a 
critical reassessment of such a possibility, recalling the elementary conditions with respect to the standard astrophysical 
background that would make it likely, showing that they are not fulfilled. Likewise, we argue that, as now well accepted, DM 
would need somewhat contrived properties to contribute significantly to the observed positron signal, even when including 
e.g. dumpiness effects. This means that most of natural DM candidates arising in particle physics beyond the standard model 
are not expected to be observed in the cosmic antimatter spectrum, unfortunately. However, this does not prevent them from 
remaining excellent DM candidates, this only points towards the crucial need of developing much more complex detection 
strategies (multimessenger, multiwavelength, multiscale searches). 
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INTRODUCTION 

Since its discovery in the early 1930'sby Zwicky [1], the DM issue has remained unsolved. There are basically two dif- 
ferent theoretical ways to address this issue, one considering a new additional component of exotic matter in the form 
of weakly interacting massive particles (WIMPs), the other involving modifications of general relativity (sometimes 
even both). Both have relevant motivations, the former from particle physics beyond the standard model [see e.g. 
for reviews] and structure formation [see a more detailed discussion in e.g. 4], the latter from more empirical attempts 
at the galaxy scale [6] or more recently in the context of extra-dimensional theories [e.g. |5fl. One of the appealing 
flavor of the former hypothesis, that we will consider in the following, is the possibility to test it with a broad variety 
of existing or coming experimental devices. Among interesting astroparticle signatures, gamma rays and antimatter 
CRs have long been considered as promising DM tracers 0Sl> but it is only recently that precision data have become 
available to look for non-standard features I19L 1 1 Oil . 

Although the rise in the local cosmic positron fraction at the GeV energy scale has been observed for a long time 
[e.g. [TJ [12, 0], the statistics recently released by the PAMELA collaboration lfl4ll is unprecedented and covers a 
much larger energy range, up to 100 GeV. The secondary origin of these positrons seems unlikely d [H, even 
if theoretical uncertainties are still large. The main questions are therefore (i) whether or not standard astrophysical 
sources may supply for such a signal and (ii) whether or not DM annihilation or decay is expected to be also observed 
in this channel. It is noteworthy that this was already discussed by Boulares lfl7ll twenty years ago, where the author 
pointed out that a pulsar origin was the best explanation to a rising positron fraction. It is not less interesting and 
sociologically striking to take a census of the articles addressing point (i) versus those focused on point (ii). Anyway, 
in this proceeding, we aim at discussing this issue concentrating on the local cosmic positron signal only, forgetting 
about other counterparts. We will first review the astrophysical backgrounds of secondary and primary origins; then, we 
will check whether DM could naturally yield prominent imprints in the cosmic positron spectrum, before concluding. 



ASTROPHYSICAL BACKGROUNDS 
Bases of CR propagation 



The global understanding of Galactic CRs at the GeV-TeV scale is rather well established. CRs are accelerated by 
shock waves at the vicinity of violent events like supernova (SN) explosions, and further diffuse erratically in the 



interstellar medium (ISM) by bouncing on moving magnetic turbulences. This diffusive motion is accompanied by 
other processes whose respective impacts depend on the cosmic ray species: convection that drives CRs away from 
the Galactic plane (negligible above a few GeV), energy losses (affecting mostly leptons), diffusive reacceleration 
(negligible above a few GeV), spallation reactions with the ISM gas (for nuclei only). The general formalism of CR 
transport was designed a long time ago in the seminal book of Ginzburg and Syrovatskii [18], and refined many times 
since then [see e.g. Il9|,|20j]. The master equation that describes the CR transport in phase-space looks like a classical 
current conservation equation: 

%J» + ® E J 1i =M. (1) 

Given a CR differential number density jY = dn/dE, the spacetime-like current is reminiscent from the Fick law 
and the heat equation, J*^ — (^V, VVK^(E) — V c \ jY}, for which the associated transport operator reads ^ = 

(d t + r spa i + Tdec, — V). The energy component is merely = J/, on which acts S>e = ds{dE/dt + Kgd]?). 
Appearing above, J2 is the source term, Kj the spatial diffusion coefficient, Ke the reacceleration coefficient, r spal / dec 

the spallation/decay rate, V c the convection velocity, and dE /dt the energy loss term. Each of these ingredients is 
by itself subject of intense researches, so that many simplifying assumptions are usually made in phenomenological 
analyses. In general, one assumes that spatial diffusion proceeds isotropically and that the diffusion coefficient is 
homogeneous in the diffusion zone, only scaling with the CR rigidity « & s . Apart from the energy losses and the 
spallation or decay rate which can be predicted independently, the propagation parameters are usually constrained 
with measurements of CR nuclei, more precisely with secondary to primary ratios like B/C [see e.g-HHIzIl]- Important 
features of such a modeling are the spatial extent of the diffusion zone (usually taken as a cylindrical slab) and the 
diffusion coefficient, and we stress that the related uncertainties are still rather large l22ll . In the following, we will 
adopt a thick cylindrical diffusion zone of radius R = 20 kpc and half-width L — 4 kpc, unless specified otherwise. 

In some cases, analytical solutions to the diffusion equation can be found in terms of Green functions Sf , which obey 
&H,E Sf^ ,£ (x <— x s ) = 5 3 (x s — x)8(t s — t)S (E s —E). For instance, assuming both steady state, which is relevant for a 
constant CR injection rate, and an infinite 3D diffusion space, thereby neglecting the spatial boundary conditions, the 
propagators for protons (or antiprotons) and electrons (or positrons) are simply given by: 

%{x^x s ) = ztt, , %{E,x^E s ,x s ) = 1 , /? exp(- ^ ^ 1 . (2) 
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Only spatial diffusion has been considered for protons (no energy losses, convection, nor spallation — fair approxi- 
mation above a few GeV). For electrons, the mixed impact of (local) energy losses dE /dt — —b(E) and of diffusion 
is encoded in the propagation scale 

l(E,E s ) = 4 / S dE'K d (E')/b(E') \ , (3) 



and other processes are neglected. This illustrates that while proton and electron transport is described by the same 
equation, these species can actually diffuse quite differently because of the relative differences in the processes they 
experience. For electrons in the GeV-TeV range, energy losses are dominated by inverse Compton (IC) scattering on 
the interstellar radiation field (ISRF) including the cosmic microwave background (CMB), and by synchrotron losses 
on the Galactic magnetic field. In the local environment, the typical energy loss timescale at Eq — 1 GeV is T ~ 315 
Myr. With a typical diffusion coefficient of Kq ~ 0.01kpc 2 /Myr, one finds X ~ 3.5 kpc, which justifies a posteriori 
the use of the local ISRF and magnetic field to compute the energy losses [1611. In the Thomson approximation, 
b{E) = (Eo/x){e = E/Eq} 2 [23], which implies that the propagation scale strongly decreases with energy. 

The previous Green functions reveal an important difference between stable nuclei and electrons: the former have 
a long range propagation scale (above a few GeV) only limited by the finite spatial extent of the diffusion zone, 
while the latter have a short range propagation scale limited by energy losses. Therefore, spatial fluctuations of the CR 
injection rate will be less important for stable nuclei (except below a few GeV, when spallation and convection become 
important) than for electrons. This implies a more local origin of high energy CR electrons and positrons, and means 
that time fluctuations in their local injection rate induces strong local effects. 



Astrophysical positrons 



Positrons of astrophysical origin can be secondaries or primaries. Secondaries are produced from spallation reactions 
of CR nuclei (mostly protons and a) with the ISM gas (H and He). Primaries can be directly produced in the intense 
magnetic field hosted by sources like pulsars [24] and further accelerated in the surrounding shocked medium, or could 
also be secondaries created from spallation processes within acceleration sites like supernova remnants (SNRs) 1 25] . 
Thus, the secondary positron source term depends on the spatial distribution of CR nuclei and of the ISM gas: 

J2 s (E,x) = 4 ^l/ dE ' d0il % ,E) B >® ' (4) 

where i flags the CR species of flux (j) and j the ISM gas species of density n, the latter being concentrated within the 
thin Galactic disk. d<3ij(E' ,E) is the inclusive cross section for a CR-atom interaction to produce a positron of energy 
E. This differential cross section w 9{E' —E)9(E' —E^Ojj/E with Eth ~ 5 GeV 01611 . so that if <p is described by a 
power law 0o£*, then £?j J (E,x) PtAn$i(E,x) Oijiij(x) / (y s — 1). Moreover, since nuclei have along range propagation 
scale and since the ISM gas does not exhibit strong spatial gradient over the kpc scale in the Galactic disk 12611 . one 
can further approximate (j) and n to their local values to get rough estimates of £} s . For instance, by considering only 
the proton-hydrogen interaction, and taking (j> p (E) = 1.3cm _2 s _1 sr _1 ((E — m p )/lGeV)~ 2 - 72 12711 . = lcm~ 3 and 
O pp = 10 mb, we end up with J2 S (E) sa 2 x 10 _27 e _2 - 72 cm _3 GeV _1 s _1 , which is actually quite close to the accurate 
calculation [163. Since the positron horizon is limited to a few kpc, the source may be modeled by an injection rate 
homogeneously distributed in an infinite thin Galactic disk of half-height h ~ 100 pc like £l s (E,5c) = 2h8(z) Qos e ^ s - 
This allows to infer the local flux inside that disk, relevant for an observer on Earth: 

j3c f „. /' , 3 . , „ * smf „ ^ pc2hQ Qs fdE s ^ rs _ cQ 0s h^ce-^ s+a -^ 



M E) = ^JdE s J i>M(E,^E.AWE,,i.) = * 2,3^, 1} • (5) 

Here, we have taken a diffusion coefficient K c i(E) = KqE 8 and an energy loss rate b(E) = (Eq/t) £ a , where a = 2 in 
the Thomson approximation. Using 5 = 0.7, a = 2 and the values given above for the other parameters, Eq. (0 gives 
<j) Si ®(E) ~ 8.5 x 10 _3 e _3 57 cm _2 GeV _1 s _1 sr _1 , overshooting the exact calculation by a factor of 2 only [ 16, 28]. 

Eq. © makes explicit the influence of the main propagationparameters: the energy loss timescale x and the diffusion 
coefficient normalization Kq set the positron flux amplitude, and their energy dependence slightly shapes the spectrum. 
Since the B/C ratio constrains mostly the ratio K(E)/L 1I22I1 . where L is the vertical extent of the diffusion zone, it is 
not surprising that the min (max) model of Donato et al. 12911 . which was designed to minimize (maximize) the primary 
antiproton flux coming from DM annihilation, is actually found to maximize (minimize) the secondary positron flux. 
Indeed, it is associated with a small value of L — 1 kpc (15 kpc), which has a corresponding small (large) value of 
Kq to fulfill the B/C constraint. Likewise, the logarithmic slope of the diffusion coefficient 8 is larger in the min 
setup, leading to a softer spectrum than in the max case. These extreme configurations, all compatible with the B/C 
constraints, are useful to bracket the theoretical uncertainties. 

In Fig. Q] we plot the latest results obtained by Delahaye et al. i28tl for the secondary positron flux, where a full 
relativistic treatment of the energy losses, i.e. beyond the Thomson approximation, was used at variance with the 
earlier calculations by Moskalenko and Strong [15] and Delahaye et al. 1 1 6j| . The right panel shows the secondary 
flux at the Earth, where the top of atmosphere (TOA) signal is corrected with a Fisk potential of 600 MV to account 
for solar modulation, while the right panel is the corresponding positron fraction defined by / = <j) e + /((j> e + +(j> e -)- 
For the fraction, we fitted the electron flux on the AMS data [30] below 20 GeV, and the full denominator itself 
on the Fermi data ll3l1l above. Predictions are shown against the data from [f32L l33i l30ll for the positron flux, and 
from IU2T, I34 13511 for the positron fraction. For the latter, are also reported the recent results obtained by the PAMELA 
collaboration 11411 . We stress that though the theoretical uncertainties are large, of about one order of magnitude 
in terms of flux, our predictions encompass the data. However, from the spectral trend observed in the positron 
fraction, it seems unlikely that the excess observed by PAMELA is of secondary origin. This naturally leads to 
the question of whether or not standard astrophysical sources may provide enough primary positrons to explain 
this fraction rise, which, from Fig. Q] should amount up to ~ 5 — 10 times the secondary flux around 100 GeV 
(0 s (lOOGeV)~4x KT^cm^GeV-'s-^r- 1 ). 

For primary electrons and positrons, the approximations made above hold, except that the source term will differ. 
Assuming again that standard sources (SNRs and pulsars) are homogeneously distributed in a thin disk of volume Vd = 
2nhR 2 , the injection of CR can be written as follows: J2 p (E,x) « 2h8(z) ^) p e^ p e^ E l Ec , where we have introduced 




FIGURE 1. Left: secondary positron flux and theoretical uncertainties. Right: corresponding positron fraction. These plots are 
adapted from Delahaye et al. 1 28] . 



an energy cut-off E c , and where the normalization =2op> which carries the dimensions, can be fixed from energetics, 
e.g. by requiring that the total rate of injected energy is set by the supernova (SN) explosion rate times the energy 
input associated with pulsars or SNRs. For SNRs, we can impose that £?op — r s „(/ E sn ) j {Vj J dE E e~Y p e~ E l Ec }, 
where F m is the SN explosion rate in the Galaxy and E m is the explosion kinetic energy of a single object whose 
a fraction of / is transferred to electrons; for pulsars, we would use instead Sq, the magnetic energy of which a 
fraction / would be converted into electron-positron pairs. Thus, the flux of primary electrons and positrons can 
also be approximated with Eq. ©, replacing for the normalization and keeping in mind that the spectral index y p 
is different, which again turns out to be a fair approximation [28]. The ratio of primaries to secondaries, given by 
R(E) = (j)p^Q (£) j (j> s ^Q (E) « (Js?op/=2?Op) £ allows to perform a quick estimate of the pulsar contribution. With 
reasonable values {y p ,f$Q,T) — (2.2, 10 47 erg,5/cy), we have /?(100GeV) > 1, which is sufficient to explain the 
positron fraction data. Anyway, even when treated more accurately such a modeling suffers much larger theoretical 
uncertainties than secondaries. First, the fraction of accelerated leptons and the averaged injected spectrum are not yet 
very well constrained by dynamical studies of sources, while important numerical efforts have been undertaken on this 
topic for a few years [e.g. 36]. Second and more dramatic, since the explosion rate of supernova? is only of a few per 
century in the whole Galaxy, the time and related spatial fluctuations become sizable locally and makes it difficult to 
justify a smooth injection rate, at least at the kpc scale around an observer. This is particularly relevant for the high 
energy component of the spectrum for which the typical propagation scale is short, and of which local sources are 
therefore expected to provide the main part. This is actually well known for decades 13711 . and was well illustrated by 
e.g. Kobayashi et al. 113811 for electrons. An important consequence of these local fluctuations is that features in the 
local spectra of CR electrons and positrons are expected. Anyway, despite the very large uncertainties, local pulsars, 
which are observed in number in the solar vicinity and whose properties can be constrained, can inject an amount of 
pos itrons that is sufficient to explain a rising positron fraction. This was recently nicely discussed in Malyshev et al. 
1 3911 . A more detailed study of primary electrons and positrons including local sources will be found in Delahaye et al. 
I 28tl . where it is shown that all current observations can be rather well reproduced with reasonable parameters. 

To conclude this part, we stress that the background to consider when looking for exotic signatures in the positron 
(or electron) spectrum is not only made of secondaries, but also of astrophysical primaries. Moreover, despite the 
large theoretical uncertainties affecting current predictions, standard sources seem capable to yield the necessary 
amount of positrons that may explain the positron fraction fairly naturally, without any over-tuning of the parameters. 
Likewise, we emphasize that the time and spatial fluctuations of the local injection rate — local sources — can lead to 
a broad diversity of features in the measured spectrum lr28Tl . which makes it difficult to disentangle different primary 
components. Finally, it seems now clear that we are far from a standard model of Galactic CRs, especially in the 
lepton channel, and many issues remain to be addressed in the future, from the CR source description to a more 
refined propagation modeling. 



DARK MATTER AND POSITRONS 



Positrons were long thought to be good tracers of DM annihilation precisely because they were expected to be of 
secondary origin only, i.e. with a low level and predictable astrophysical background. As argued in the previous 
section, this statement is likely not valid anymore. Anyway, to keep the reasoning as general as possible, let us recall 
some basic conditions for a cosmic messenger to be a good tracer for any exotic signal: (i) the background is not too 
high with respect to the expected signal, given an experimental sensitivity; (ii) the background is known or predictable, 
and controled; (iii) specific spectral features in the signal make it unambiguously distinguishable from the background. 
We did not yet discuss condition (i), but it is clear from the previous section that conditions (ii) and (iii) cannot be 
fulfilled. From this simple argument, we can hardly hope, at least with current data, to identify a clean DM signature 
in the local positron spectrum. Nevertheless, it is still well-grounded to ask whether DM is about to provide a sizable 
signal, should it be mixed with other components. If so, there might still be some hopes for isolating it with future 
experiments, provided the astrophysical background is being much better understood in the meantime. 

We will first review the predictions that can be derived in general cases when modeling the Galactic DM halo with 
a smooth distribution. Then, we will discuss the potential impact of DM substructures on the expected signal. We will 
focus our discussion on annihilating DM, disregarding decaying candidates. 



The smooth approximation 

Structure formation in a ACDM universe involves non-linear processes as soon as the linear growth of perturbations 
triggers the gravitational collapse of objects. In this theoretical framework, galaxies are expected to have formed 
around redshift z ~ 6, and have consequently left the linear regime for a long time, so that only numerical experiments 
can provide detailed information on the DM distribution in those objects. The advent of high resolution numerical 
simulations in the study of structure formation has led to major breakthroughs during the last two decades, providing 
a fairly good understanding of the properties of the large scale structures that are observed in current surveys. There 
are still some important mismatches at the Galactic scale, but this might be due to the important impact of baryons 
which have not been included to those simulation until recently [see e.g.|40j for a discussion on the small scale issues]. 
Anyway, disregarding the potentially large effect induced by baryons, it seems that DM structures at the galactic scale 
are predicted to have similar smooth DM density distributions, almost spherical and scale-independent. Such a generic 
density profile was named NFW after their authors Navarro, Frenk and White l4lll . and is typified by the following 
function: 

p(0=P.v n ^ / / r \ 2 , (6) 
(1 +r/r s y 

where r s is a scale radius beyond which the logarithmic slope goes from -1 to -3; p s is the scale density. It is noteworthy 
that more recent results give quite similar smooth DM components scaling like p ~ (r/r s )^ in the central regions, with 
y ~ 1, though the lack of resolution prevents from making clear predictions at the very center [see e.g. H2,|43l]- Using 
such a density shape, derived from theoretical constraints, to describe our Galaxy needs to account for additional 
observational constraints from stellar kinematics. Basically, although the baryon modeling comes into play with 
uncertainties, a local density of Pq(Rq = 8 kpc) ~ 0.3GeV.cm~ 3 associated with a scale radius of r s < 20 kpc are 
fairly compatible with current data [e.g. |44|, |45l |46[] . Now, since the injection rate of the DM annihilation products is 
« p 2 , the precise values of the logarithmic slope may have a strong impact on predictions, especially for annihilation 
close the Galactic center (GC). This is particularly true for 7-ray flux predictions 14711 . but not really for CR positrons 
in the GeV-TeV energy range, since, as discussed above, they cannot pervade beyond a few kpc — for longer range 
antiprotons, the signal coming from there is anyway diluted by diffusion. Uncertainties in the flux amplitude will be 
therefore mostly set by those on the local DM environment. 

To summarize, let us write the source term associated with DM annihilation: 

I P© ) dE 2 {m x j 



where 5 = 1/2 (1) for Dirac (Majorana) fermionic WIMPs, and 1 for scalar WIMPs; (<7v) is the WIMP annihilation 
cross section, m x the WIMP mass and dN(E) the number of CR positrons injected in the energy range dE. The 




FIGURE 2. Letf: Predictions of the primary positron fluxes for a smooth NFW DM halo and 3 different annihilation final states, 
as compared with the secondary background. Right: Boost factors obtained by Pieri et al. 15111 for direct annihilation in e + e~ as 
calculated in the frames of the Via Lactea II 1 43] and Aquarius li^l setups. 



different WIMP candidates should have similar annihilation cross sections of (<7v) ~3x 10~ 26 cm 3 /s if they decouple 
thermally from the primordial bath as a consequence of expansion in the early universe, such a value being fixed by the 
present cosmological DM density [the generic method to compute the relic abundance can be found in 48]. The local 
positron spectral shape will primarily depend on the annihilation final states. Three typical final states may basically 
typify the main features of the positron injection spectra associated with DM annihilation or decay, (i) quarks, say bb, 

(ii) and (iii) l + l~ with any charged lepton, say e + e~ . The spectrum is getting harder and harder from (i) to 

(iii) . 

To check whether usual WIMP candidates are about to give a sizable positron flux, it turns out useful to derive the 
flux in the asymptotic limit of very short propagation scale, which is a rough approximation valid at high energy. The 

corresponding propagator is therefore @ e (E,x «— E s ,x s ) ► S(E S — E)8 3 (x s —x)/b{E) — which only differs from 

the diffusionless limit by the term 8(E S — E) that ensures X — » here. If we further assume a direct annihilation in 
e + e~, so that dN(E)/dE = 8(E — m x ), then the asymptotic and exact flux limit reads: 
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where we have used the Thomson approximation for the energy losses. Note that for m x = 100 GeV, this result is 
pretty close to the prediction of the secondary positron flux at 100 GeV. Since the positron fraction measurement at 
100 GeV implies ~ 5 — 10 times more primaries than secondaries, this means that boosting the local DM density by a 
factor of ~ 3 is enough to feed the PAMELA data rather significantly. Nevertheless, this is, at least to our knowledge, 
the unique example for which one may recover the observed positron fraction at ~ 100 GeV without over-tuning the 
annihilation cross section. Indeed, for other annihilation final states and/or larger WIMP masses, one needs to boost 
the signal by 2 to 4 orders of magnitude to get enough positrons at 100 GeV to match the measurements fg. g. I49L 15011 . 
This is illustrated in the left panel of Fig. [2] where we compare predictions assuming the three annihilation final 
states discussed above, a smooth NFW density profile and the parameters used in Eq. ([8]). In this plot, the secondary 
background is taken from Delahaye et al. 1 16], derived in the Thomson approximation for the energy losses. 

To summarize this basic analysis, it seems that few WIMP candidates with thermal relic abundance may provide 
sufficiently positrons to feed the positron fraction data naturally. Indeed, only those WIMPs with masses around ~ 100 
GeV and with direct annihilation in e + e~ do not need arbitrarily high boost factors with respect to a smooth description 



of the density profile. Needless to say that there are poor motivations for such models in particle physics beyond the 
standard model — couplings to heavier leptons would lower the positron yield, however, why 100 GeV particles 
should couple only to e + e~ ? — and that there might already exist limits coming from colliders. Therefore, it is fair 
to conclude that DM annihilation does not provide a natural explanation to the PAMELA data. However, it is not less 
fair to ask about the potential impact of DM substructures that are predicted in the frame of ACDM, and that could 
enhance the local DM density. Note that we have studied the impact of relaxing spherical symmetry for the smooth 
halo profile in Lavalle et al. 15211 . showing that this has poor effect on predictions. 



Clumpiness effects 

DM substructures (called also subhalos or clumps) are predicted in the ACDM paradigm and observed in cosmo- 
logical N-body simulations. The smallest scales that can grow and further collapse are those encompassed within 
the WIMP free streaming scale which is set by their intrinsic properties (mass, couplings). For generic WIMPs, the 
minimal mass scale ranges within lO" 11 - 1O~ 3 M [53]. The mass function of these subhalos is usually found close 
to the oc M~ 2 prediction of the Press-Schechter theory of self-similar gravitational collapse l54l l55ll in cosmological 
N-body simulations [for recent results, see l43l l42tl . Consequently, without loss of generality, the subhalo distribution 
in a Milky-Way-like object may be written as B56I1 : 

dn c \(m c \,x) _ N ^ d^mjm^x) d^ v (x) 
dm c \ dV dm c \ dV 

where N c \ is the total number of subhalos, and where d£? m and denote the mass and spatial probability distribution 
functions, which are normalized to unity over the whole galaxy extent. N c \ can be constrained from N-body simulation 
results, at least in the available resolved mass range — the most recent simulations involving billions of particles can 
resolve clumps down to 10 3 — 10 4 M Q at the galaxy scale 11431 14211 . When extrapolated down to a minimal scale of 
1O~ 6 M for subhalos, this number is found in the range 10 14 — 10 16 in a Milky-Way-like galaxy ll56l ETll . Analytical 
studies on the tidal disruption of these subhalos in galaxies due to gravitational interactions with the disk or stars 
show that an important fraction may actually survive [57]. Besides, although rather spherical, the spatial distribution 
of subhalos turns out to be different from the smooth DM profile p (r) since it exhibits a core radius rather than a cusp 
in the central parts of N-body galaxies, which is likely due to efficient tidal disruption, and a r~ 2 form rather than 
r~ 3 in the outskirts. Such a behavior is called antibiased, because (d£?v/dV)/p(r) r HH], but it is still not clear 
whether it is still valid for lighter objects, far from being resolved in N-body simulations. Finally, the mass distribution 
should in principle depend on the location in the galaxy to account for tidal disruption. One can include this effect by 
calculating the maximal subhalo mass as a function of the galactic radius, keeping constant the normalization of the 
mass function, which is performed in the full mass range 05111 . Anyway, as mentioned above, dS? m °< m c \~ a dm c \ with 

Since substructures are generically predicted in the ACDM paradigm, and are expected to be impressively numerous 
in galaxies if DM is made of WIMPs, it is important to include them for consistent calculations of astrophysical signals. 
Indeed, since the DM annihilation rate is proportional to the squared DM density, the presence of subhalos in the local 
environment can have strong impact on the antimatter flux (as well as on the diffuse photon emission). We have 
therefore to derive a method to add the subhalo contribution to the smooth one. Since we have sketched a phase-space 
distribution of subhalos in Eq. (O, we may think about treating the flux coming from a single object like a stochastic 
variable, which actually turns out to be correct and powerful 11591,15611 : the typical range of subhalo scale radii, where 
most of the annihilation proceeds, is indeed smaller than the typical propagation scale, so that clumps can be safely 
treated as point-like sources. We still further need to specify the properties of a single object, i.e. its mass m, position 
R in the Galaxy, inner density profile and the amount and spectral shape of positrons that it injects. It is conventional 
to define the subhalo extent by the radius r v at which the average subhalo density is 200 times the critical density 
p c of the universe today. Even when fixing the shape of the inner profile, taking e.g. the NFW model, constraining 
the associated scale parameters r s and p s is in principle much more complicated, since they depend on the formation 
history. Nevertheless, this history shows up an evolving correlation between the concentration, defined by c v = r v /r s , 
and the subhalo mass — the less massive the more concentrated because formed earlier, in a denser universe. Thus, 
the knowledge of this concentration function at z = allows to specify the subhalo parameters entirely. Not only 
does this concentration function depend on the subhalo mass, but also on the its location in the Galaxy, since more 
concentrated objects resist more efficiently to tidal effects. It is convenient to define the annihilation volume of a single 



object: ${m,R) = 4x$drt 2 { Pd ^ r) } . This actually defines the volume needed from a constant density of p 
to produce the actual subhalo injection rate, and somehow measures the ratio of its intrinsic emissivity to the local 
emissivity. Armed with this definition, it is straightforward to derive the local average flux associated with the whole 
subhalo population l59l l56ll : 

c i©( £ ) = ^N d y{{t;{x s )),J(x Q ^x s )) v where &{x Q ^x s ) = j™* dE^^E ^x s ,E s ) dN[E] 
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where Q x denotes the average performed with the distribution and the latest line is the limit corresponding to a 
vanishingly small propagation scale A — > 0. 

In order to check whether a clumpy DM halo leads to a larger positron flux, it is useful to compute the ratio of 
both predictions. Of course, clumps are not merely additional mass in the halo, there should be some consistency 
as well as observational constraints to obey. In fact, the census of subhalos in N-body simulation rests on the mass 
resolution: since the galactic scale is described quite accurately, one can consider that whatever the discreteness of 
its content, an N-body galaxy will have a constant mass. Therefore, to model the Galaxy in a consistent manner, a 
certain fraction of mass should be removed from the smooth DM profile when adding clumps: p(r) — > (1 — f)p(r), 
where / = N c \(m c \) m /My^. The potential problem with such a procedure is that as soon as the spatial distribution 
of subhalos differs from the smooth distribution, then the average local DM density is modified [56], which can lead 
to comparing situations with different local average DM densities. Anyway, we can now derive the ratio of the flux 
prediction for the smooth case to that for the clumpy case, so-called boost factor. 

0© (E) dV 

where the limit of vanishingly small propagation scale is obtained from Eq. ([8]l and Eq. ( fTOb . This expression is quite 
natural: the boost limit is only given by the local number density of objects N c id^'y(xQ)/dV times the average 
annihilation volume of a single object (£,(x & )) m (which is normalized to the local smooth luminosity by definition, 
see above). Before taking a numerical example, we emphasize that the general expression of the boost factor is a 
function of energy. Indeed, for large propagation scales, i.e. low energy, the signal coming from the cuspy smooth 
distribution in the GC will certainly dominate the total flux, while subhalos may dominate at short propagation scales. 
This is an important feature which is very often neglected. Notice also that the global subhalo flux is associated with a 
statistical variance that increases as the number of objects decreases in the relevant propagation volume: this variance 
does therefore increase with energy l59ll56ll . Correspondingly, the boost factor can also have a large variance provided 
subhalos dominate over the smooth contribution, otherwise it is diluted. 

Let us take a very simple and rather optimistic example in which we assume a total number of iV c i = 10 17 subhalos 
with fixed masses 1(T 6 M Q , with inner NFW profiles, a fixed concentration of c v = 100, and distributed according 
to a cored isothermal profile of core radius r c = 20 kpc, extended up to 280 kpc in the Galaxy taken with a mass of 
1O 12 M . We have therefore d^ v (r@ = 8kpc) /dV «5x 10^ 7 kpc" 3 , a subhalo mass fraction of / = 0. 1, so that: 



Stf « (1 -ff +Na 200^(3 + C; (3 +Cv )) = 

This result is quite modest even with rather optimistic parameters {£, °< elm). Should have we taken a subhalo spatial 
distribution tracking the smooth component, we would have found 4. From this simple estimate, it seems 

unlikely that light subhalos, even if numerous, provide a strong enhancement of the positron signal in average. 
Interestingly, this naive reasoning gives a result which is actually quite close to accurate calculations involving more 
complete subhalo models 1561 , 15 ill . For illustration, we show in the right panel of Fig. [2]the boost factors and associated 
statistical variances obtained from the subhalo settings derived in Diemand et al. B43I1 and Springel et al. B42I1 . where 
WIMPs of 100 GeV annihilating in e + e~ are assumed. It is noteworthy that in this plot, boost factors are shown to 
be in fact less than 1 ! This is a consequence of the decrease in the local density due to the addition of subhalos, 
constrained so to keep the Galaxy mass constant. The statistical variance shown as colored bands expresses the fact 
that few nearby objects can actually dominate the whole subhalo contribution. This would induce peculiar features in 
the positron spectrum, hardly distinguishable from those predicted for nearby standard astrophysical sources [e.g. 60]. 



CONCLUSION 



We have argued that the rising local positron fraction observed by the PAMELA satellite is unlikely of secondary 
origin. Then, we have discussed the potential yield from astrophysical sources of primary positrons, emphasizing 
local pulsars as our best candidates, and stressing that those conclusions were already sketched twenty years ago 11711 . 
Although the overall electron and positron data can be very well explained with standard astrophysical processes, 
we have finally stressed that we are still far from a standard model of cosmic rays, since the current theoretical 
uncertainties on propagation, source modeling as well as ISM modeling make it rather unfair to claim for clean 
predictions; at the moment, we can only provide rough estimates still tuned — though reasonably — to reproduce 
the data. In light of this discussion, it seems to us that the rising positron fraction, as well as the so-called electron 
excess sometimes seen in the Fermi data, are no longer theoretical issues, since standard and not contrived explanations 
are available. Remains wide open the question of identifying and modeling accurately the local sources of cosmic ray 
leptons in order to sustain this solution on more detailed grounds. These are rather good news for this research domain, 
which, besides the need for important theoretical efforts, can benefit a copious amount of experimental data. Of course, 
this implies a multimessenger and multiwavelength analysis, which are mandatory for consistency purposes. 

Regarding the DM hypothesis, we have shown that usual WIMP candidates are not expected to contribute sig- 
nificantly to the local positron flux, even when treated in a self-consistent framework including subhalos. The only 
possibility without over-tuning the annihilation cross section allowed for thermal relics is to consider direct produc- 
tion of e + e~ and a mass scale of ~ 100 GeV, which is not motivated in particle physics theories beyond the standard 
model. Likewise, we have also stressed that should DM yield a sizable positron signal, it would be difficult to disentan- 
gle it from standard astrophysical sources. This illustrates the fact that the basic conditions that would make positrons 
good DM tracers are not fulfilled: not only is the background larger than the signal, but, more important, it is not yet 
under control. 

Anyway, we underline that though unlikely contributing to the local positron flux, WIMPs remain excellent DM 
candidates. The crucial issue of their detection is still challenging, since their expected properties have made them 
continuously escape from observation despite the advent of important experimental devices, especially in high energy 
astrophysics. It seems important to develop more complex strategies based on multi-messenger, multi-wavelength 
and multi-scale approaches, in which large efforts should be made to quantify and minimize the associated theoretical 
uncertainties. Other detection methods are also very important, among which the LHC results are particularly expected. 

Acknowledgments: It is a pleasure to thank T. Delahaye and R. Lineros for sharing a great deal of work in this 
topic. 
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